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ABSTRACT 



Previous comparisons of optical and X-ray observations of clusters of galaxies have 
led to the so-called - discrepancy" that has persisted for the last decade. The standard 
hydrostatic-isothermal model for clusters predicts that the parameter Papec = (^r / i^'^ / t^''^p) M 
which describes the ratio of energy per unit mass in galaxies to that in the gas, should 
equal the parameter (3fit (where Pgasir) oc Pgaiif)^^^^^) determined from the X-ray surface 
brightness distribution. The observations suggest an apparent discrepancy : [3 spec ~ 1.2 
(i.e., the galaxies are "hotter" that the gas) while (3 fa ~ 0.65 (i.e., the gas is "hotter" and 
more extended than the galaxies). Here we show that the discrepancy is resolved when 
the actual observed galaxy distribution in clusters is used, Pgaii^) oc 7--2-4±o.2^ instead 
of the previously assumed steeper King approximation, Pgai{r) oc r~^. Using the correct 
galaxy profile in clusters, we show that the standard hydrostatic-isothermal model predicts 
Pspec = Pfit — (1.25 ± 0.10) X (3fit, rather than Pspec — Pfit (where Pfu is the standard 
parameter using the King approximation, and is the corrected parameter using the 
proper galaxy distribution). Using a large sample of clusters, we find best-fit mean values 
of Pspec = 0.94 ± 0.08 and (3^^^ = 1.25 x (Sfu = 0.84 ± 0.10. These results resolve the /3 - 
discrepancy and provide additional support for the hydrostatic cluster model. 

Subject headings: Galaxies: clustering - X-rays: galaxies 
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1. INTRODUCTION 



The gravitaional mass of clusters of galaxies has traditionally been estimated from 
the dynamics of the cluster galaxies using the virial theorem, revealing large amounts of 
dark matter (Zwicky 1933, Peebles 1980). More recently. X-ray emission from the hot 
intracluster gas has been used to estimate cluster masses by utilizing the intracluster gas 
temperature and the density profile as a tracer of the cluster potential assuming hydrostatic 
equilibrium for the cluster gas (Cavaliere & Fusco-Femiano 1976, Bahcall & Sarazin 1977, 
Forman & Jones 1984, Sarazin 1986, Hughes 1989, Evrard 1990, BahcaU & Cen 1993). 
The two methods yield comparable masses. However, one problem has persisted over the 
years : the so-called - discrepancy" of clusters. This problem reflects the discrepancy 
between the observed parameter /3spec , determined from the X-ray temperature and cluster 
velocity dispersion, and Pfu , determined from the gas versus galaxy density profiles in 
the clusters. If the standard hydrostatic cluster model is correct, the two parameters are 
expected to have similar values; however, observations suggest that Pspec ~ (1.5 — 2) x 
(Sarazin 1986, Evrard 1990). This discrepancy has been an unsolved puzzle in the study 
of clusters of galaxies for nearly a decade. 

In this paper, we show that the "/9 - discrepancy" results mainly from assuming a 
galaxy distribution in clusters which is too steep (the King approximation), thus causing 
misleading conclusions. The - discrepancy" is resolved if the actual observed profile of 
the galaxy distribution in clusters is used. 



The standard model for the structure of mass in clusters assumes that both the gas and 
the galaxies are in hydrostatic equilibrium with the binding cluster potential (Cavaliere & 
Fusco-Femiano 1976, Bahcall & Sarazin 1977, Forman & Jones 1984, Sarazin 1986, Evrard 
1990). In this model the gas distribution obeys ^ = ^ = -^^P-, where Pgas 

and Pgas are the gas pressure and density profiles, (f) is the cluster potential, and M(r) is 
the total binding cluster mass within radius r of the cluster center. The cluster mass can 
thus be represented as 



2. THE (3 - DISCREPANCY 





) r 



(1) 



dlnr dlnr 
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where T is the intracluster gas temperature, and fxrup is the mean particle mass of the gas. 
The galaxies in the cluster respond to the same gravitational field, and thus satisfy 

= _^ (dJl^El^ + ^ + 2A)r (2) 
G \ dlnr dlnr ' 

where Ur is the radial velocity dispersion of the galaxies in the cluster, Pgai{r) is the galaxy 
density as a function of r, and A represents a possible anisotropy in the galaxy velocity 
distribution [A = 1 — (at/cTr)'^, where t and r represent the tangential and radial velocity 
components] . 

Relations (1) and (2) thus yield : 

p _ o'r din p gasir) / dlnr + dlnT / dlnr 

spec I dlnpgai (r) I dlnr + dlna^ / dlnr + 2 A 

where the Pspec parameter is defined by the left side of relation (3). This spectral /? 
parameter can be determined directly from observations of cluster velocity dispersions and 
gas temperatures; it represents the ratio of energy per unit mass in the galaxies to that in 
the gas and is observed to be ~ 1 (see below). (Here and below we assume that the cluster 
velocity dispersion is isotropic, thus the radial velocity dispersion ar is comparable to the 
observed line-of-sight velocity dispersion a). The right-hand side of equation (3) relates 
the (3spec parameter to the density profiles of the gas and galaxies in the cluster as well as 
to the temperature and velocity profiles. 

For isothermal galaxy and gas distributions and isotropic galaxy velocities (a model 
that is traditionally used to describe clusters), one has dlnT/dlnr = dlna^/dlnr = ^ ~ 0; 
relation (3) then becomes 

^ ^ ^ dlnpgas{r)/dlnr _ 

spec kT/jimp din p g ai{r) / dlnr 

the f3spec parameter determined from ar and T (left side of eq. 4) should approximately 
be equal to the ratio of the slopes of the gas to galaxy density profile, defined as |3'j^^ 
(right side of equation 4). The solution of equation (4) yields Pgas{f) oc Pgai{r)^f^* and 
Pspec — Pfit ■ The parameter Pj^^ is determined by fitting the observed X-ray surface 
brightness distribution in clusters. The galaxy distribution needed to normalize the Pj^^ 



determination has been generally assumed to be the King (1962) approximation, Pgaiir) = 
Pgai{0){l + {r/RcY) , where Pgai{r) oc for large radii (i.e., r > Rc, where Rc is the 
cluster core radius). Therefore, Pgasir) = pc,as(0)(l + {r/RcY) ^""^ ; here (3fit without 
the superscript "c" for "correct" denotes the standard definition of (3 fa using the King 
approximation (see, e.g. Abramopoulos & Ku 1983, Jones & Forman 1984, Sarazin 1986, 
Evrard 1990, and references therein). Therefore, Pfn has previously been determined from 
the relation 

Pgasir) [a Pgasirff''] OC r"^^^" (5) 

where Pfn represents the determination using the King approximation, and P^^^ is the 
"correct" value from equation (4). 

X-ray observations (Abramopoulos & Ku 1983, Jones & Forman 1984, Henry et al. 
1993) yield Pfn values in the range ~ 0.5 to 0.9, with a median of Pfn — 0.67 ±0.10 (rms); 
this corresponds to a gas density profile of 

Pgasir) oc r-2°±0-3 Rc<r< 1.5 h'^ Mpc (6) 

On the other hand, observations yield Papec = cr^ / {kT / finip) values in the range of ~ 0.5 
to ~ 2 with a median Pspec 1 - 1.4 (Mushotsky 1984, Sarazin 1986, Evrard 1990, Edge & 
Stewart 1991, Lubin & Bahcall 1993). This inequality between Pspec and Pfu contradicts 
the expectation from the hydrostatic model. We discuss below how this discrepency is 
resolved. 

3. RESOLVING THE - DISCREPANCY 

The "/3 - discrepancy" is based on the assumption that the galaxy distribution in clus- 
ters follows the King approximation, Pgai{r) oc at large radii. This assumption, how- 
ever, is inaccurate. Specific studies of the galaxy distribution in clusters yield a shallower 
profile. The galaxy-cluster cross-correlation function, which represents the average net 
galaxy density profile around rich clusters, yields Pgai{r) oc r~^'^ (Lilje & Efstathiou 1988), 
or Pgai{f) oc r"^'"^ (Seldner & Peebles 1978, Peebles 1980). The average profile of a sample 
of rich clusters yields Pgai{r) oc r~^-^ (Bahcall 1977) (all for 0.5 < r < 1.5 Mpc). 
The profile becomes even shallower, Pgai{r) oc for r < 0.5 Mpc (above references; 
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see also Beers & Tonry 1986). We use below the observed range of the average galaxy- 
density distribution in rich clusters 

Pgalir) oc r-2-^±°-2 0.5 h'^ < r < 1.5 h'^ Mpc (7) 

Inserting relations (6-7) in equations (4-5) we find : 

Pgasir) OC Pgalir^f^* OC r-(2-4±0.2)/3J,, ^ ^Sfifu ^ ^-2.0±0.3 (g) 

where the last term in (8) represents the observed X-ray profiles (Sect. 2). The corrected 
Pj^^ parameter introduced here is related to the parameter determined using the King 
approximation, Pfn (eq. 5), by 

^f'' 2.4 ±0.2 ^f'' = ^^-^^ ^ ^f'' 

Since the X-ray observations yield a mean value of Pfu = 0.67 ± 0.02 (Jones & Forman 
1984, Henry et al. 1993), the corrected Pj^^ parameter (eq. 9) has a mean value of 

Pjit = (1.25 ± 0.10) X Pfit = 0.84 ± 0.10 (10) 

If the fit is dominated by X-ray emission at small separations, where Pgai oc r~^, then (3j--^ 
~ 1. 

Detailed X-ray and optical observations of the cluster A2256 (Henry et al. 1993) 
using recent X-ray data from ROSAT support the above analysis; the observations yield 
Pgas oc r~^-^, Pgai OC r"^'^, Pfu = 0.795, and therefore Pj-^ = 0.96, fully consistent with 
relations (7-10). 

A recent analysis of the largest-available sample of clusters with measured T and line- 
of-sight velocity dispersions a (each having more than twenty galaxy redshifts per cluster) 
yields a best-fit Pspec parameter for the entire sample (Lubin & Bahcall 1993) : 

2 2 

ki / iiirip ki I iiirip 

This best-fit {3 spec was obtained from a fit to the observed (t{T) relation for the sample 
[cr = {Pspec/ piT^p)^'^ (kT)^-^] using the average observed fj, = 0.58 for clusters (Edge & 
Stewart 1991). The median value for the sample is Pspec {median) = 0.98. 



The directly observed values of (3 spec and (eqs. 10-11) are consistent with each 
other as expected from the hydrostatic- isothermal model (eq. 4). Various effects such 
as possible substructure in some clusters (Geller & Beers 1982), velocity anisotropy or 
contamination, cooling flows, or incomplete thermalization of the gas may contribute to 
the observed scatter (see Figures below) but do not significantly affect the main results 
obtained above (see Lubin & Bahcall 1993). 

The agreement between figpec and persists even if the distributions are not isother- 
mal. Some clusters show a decrease in their temperature and velocity dispersion profiles 
(e.g., Coma; Kent & Gunn 1982, Jones & Forman 1992). Using relation (3), with temper- 
ature and velocity gradients of dlnT/dlnr ~ dlna^/dlnr ~ —0.5 to -1 (as appears to be 
suggested by the results in the above references), as well as allowing for a small amount 
of velocity anistropy in clusters {A < 0.2), we find f3j-^ ~ 0.9 — 1 (from the right side 
of eq. 3). This value is in excellent agreement with the directly measured value of (3spec 
= 0.94 ±0.08 (eq. 11). 

The results obtained here are illustrated in Figures 1-3. In Figure 1, we plot the 
observed Pj^^ and Pspec values for all clusters where both parameters have been measured 
(Jones & Forman 1984 for Pfn ; Lubin and Bahcall 1993 for Pspec )■ The original Pfu 
values (assuming the King approximation) are presented in Figure la, while the corrected 
Pf^^ values are presented in Figure lb. While Figure la clearly refiects the /3 - discrepancy 
(i.e. Pspec > Pfit)j this discrepancy is eliminated in Figure lb, where the corrected /^^^^ 
values (corresponding to a galaxy density profile of pgai oc 7--2-4±o.2-j ^^.^ plotted. The 
ratio Pspec / Pfit presented for all clusters in Figure 2. It is clear from the results 
in Figure 2 that, within the observational uncertainties, no (3 - discrepancy is apparent 
{Papec /Pfit ~ !)• Finally, we plot in Figure 3 the recently observed gas and galaxy density 
profiles in the cluster A2256 (Henry et al. 1993). The data supports the main conclusions 
reached above : the gas and galaxy distributions follow each other (i.e., ~ 1) within 
the observational uncertainties. Both distributions can be represented by an p oc r~^'^ 
profile. The spectroscopic (3spec parameter of the cluster is (3spec = 1-33 ± 0.24 (Lubin & 
Bahcall 1993); these values are included in Figures 1-2 and are consistent with the entire 
set presented. 
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We conclude that no significant (3 - discrepancy exists in clusters of galaxies. The 
hydrostatic model therefore provides a consistent fit to both the X-ray and the optical 
data. 

We would like to thank J. P. Ostriker, P.J.E. Peebles, D. Richstone, C.L. Sarazin, and 
D.N. Spergel for helpful discussions. 
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FIGURE CAPTIONS 

Fig. 1 

(a) The observed values of the parameters Pfu (using the King approximation; Jones & 
Forman 1984) and Pspec {— cr"^ / kT / iinip] Lubin & Bahcall 1993) for all clusters for which 
both parameters are measured. The "/3 - discrepancy" effect, i.e., (3 spec > Pfit, is apparent 

(b) Same as (a), but for the corrected values, assuming the average observed galaxy- 
distribution Pgaiir) oc (Sect. 3). No "/3 - discrepancy" is apparent. (The scatter 
in Pspec is larger than in , as expected due to larger observational uncertainties.) 

Fig. 2 

The ratio Pspec / Pfu available clusters (same data as Figure lb). No apparent 

discrepancy is seen (i.e., (3spec / Pju ~ !)• Different symbols represent different sub-samples 
: filled squares - clusters in superclusters; triangles - isolated clusters; open squares - 
clusters at |&| < 20° (see Lubin & BahcaU 1993). 

Fig. 3 

Projected galaxy density distribution (open circles) and projected gas density distribution 
(solid curves represent ±1 sigma results) for the cluster A2256 (Henry et al. 1993). The 
galaxy density is in units of galaxies per deg^; the gas density is in arbitrary units. Both 
profiles follow p(r) oc r~^-^ outside the cluster core. 
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